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ABSTRACT 






An all-thruster three-axis stabilized attitude control 
system has been designed for the Naval Postgraduate School 
(NPS) satellite bus, ORION. The satellite is a cylinder, 19 
inches in diameter, 35 inches in length, 250 pounds maximum 
mass with 32 pounds for payloads. ORION will be ejected from 
an extended Get-Away-Special (GAS) canister. Launch from any 
GAS can configured expendable booster or the space shuttle is 
assumed. The minimization techniques of Pontryagin have been 
used to derive control laws that support fuel efficent 
operation. A minimum time cost function is applied in the 
acquisition phase to reduce rates to acceptable levels. A 
weighted minimum fuel-time cost function is used during the 
on-station phase. Bang-Off-Bang control with two switching 
curves is employed outside of a boundary region. Inside the 
boundary region, four-pulse limit cycle control with time 
constants on the order of 100 seconds is applied. Lifetimes 
in the range of 4 to 28 months are obtained utilizing these 
techniques . 



iv 



D'J-r., • 

II A 'i 
MOlUn 




THESIS DISCLAIMER 



The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I . INTRODUCTION 



The current ORION design uses spin stabilization, but a 
truly general purpose satellite should provide attitude 
control design options for possible payloads. In surveys of 
likely Navy and Air Force Space Test Program (STP) payloads 
during the past decade, the majority of the requests were for 
three-axis stabilization as shown in Table 1-1 [1]. This 
bears out for the potential customers who have shown an 
interest in the Orion bus. The intent of this thesis is to 
investigate a mass expulsion three-axis stabilization design 
as an ORION bus option. 

A three-axis stabilized Orion bus will undoubtedly be 
more expensive than a spin stabilized version. Limited 
changes to the original design are desirable. Modified 
deployable solar panels will be necessary to achieve the 
required end-of-life power. Additionally, a revised sensor 
package with more complex sun, horizon sensors and rate gyros 
will be reviewed. This thesis investigates minimum fuel-time 
optimization with respect to thruster operation using 
Pontryagin's Minimization Principle. The effects of external 
disturbances including gravity gradient and aerodynamic drag 
cannot be ignored. Models of a three-axis ORION were 
simulated. 
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Table 1-1 SUMMARY OF TYPICAL NAVY/STP 
REQUIREMENTS [1] 



PAYLOAD 



Mass 


32 lbm 


Volume 


2.36 ft 3 


Power 


34 watts 


DataRate 


5000 bits/sec 


Orbit 


200-800 nm circular 


Inclination 


0 ° -30 ° or 60 0 -120 0 


Instruments 


Particle counter or Lens 


Control 


3 Axis ,+.75' 
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II. BACKGROUND 



Today most satellites built in Western countries are 
large, complex and multi-mission/payload capable. The cost 
of additional payloads has proven exponential due to many 
factors, of which the most debilitating is that contractors 
are not motivated to organize their program structures 
efficiently. Other factors include: expensive launch from 
only two available sites, the loss incurred when a booster or 
satellite subsystem fails and the increased testing required 
before launch. This puts the free world in a highly 
undesirable position. Reserve assets and the ability to 
produce large numbers quickly i.e., surge capability for these 
high value satellites is extremely limited, if not non- 
existent. In addition, most programs have been stalled at the 
1970-75 technology level due to a general fear of using 
unproven components when the probability of success is already 
questionable and the cost is unacceptable. The concept of 
small, cheap satellites should not work against but with 
current programs. A multi-site launch capability exists from 
ships, submarines, aircraft and mobile platforms. New 
technology can be implemented and space qualified on low value 
vehicles. The ability to build up reserves and to launch 
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several satellites at a time becomes viable. The advantages 
are myriad. 

A. CURRENT SMALL SATELLITES 

The present U.S. menu of small proven space platforms 
consists of SPARTAN, NUSAT, GLOMR and LDEF [1]. SPARTAN is 
a retrievable rectangular platform launched from the Space 
Shuttle that can fly alongside for only 2.5 days with a 
payload limit of 1000 pounds. NUSAT (Northern Utah Satellite) 
is also Shuttle deployed but from a Get-Away-Special (GAS) 
canister with a specially designed launch mechanism displayed 
in Figure 2-1. 




Figure 2-1 Get-Away-Special Canister [1] 
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GLOMR (Global Message Relay) was a DARPA (Defense Advanced 
Research Project Agency) payload which used the NUSAT design. 
LDEF (Long Duration Exposure Facility) is basically a large 
(up to 20,000 pounds) , multi-experiment platform that provides 
a facility for extended space exposure. Figures 2-2, 2-3 and 
2-4 illustrate SPARTAN, NUSAT and GLOMR. Though these 
represent advancements in small satellite development, they 
are fundamentally free-flying structures incapable of transfer 
to higher orbits or long duration attitude control [1]. 
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Figure 2-3 Weber State College Second Generation NUSAT [1] 
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Figure 2-4 DARPA GLOMR Satellite [1] 
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The ability to produce small multi-purpose satellites has 
existed since the dawn of the space age. The need to produce 
a simple yet reliable, low cost satellite bus that is 
available to commercial and military interests alike has never 
been as great as it is now. The Naval Postgraduate School 
(NPS) satellite bus concept not only encompasses propulsion 
and attitude control subsystems, but also standard electric 
power, telemetry, thermal and computer subsystems. 

B. ORION HISTORY 

The history of NPS satellite development began in spring 
1985 when the then Chairman of the NPS Space Systems Academic 
Group, Dr. Allen Fuhs, taught the annual spacecraft design 
course. The project for that particular class of space 
engineering and operations students was the design of a small 
submarine launched satellite [1]. The course was taught two 
more years before Dr. Fuhs retired from NPS. Each of these 
two classes worked on applications of the satellite. One 
student, LCDR Austin Boyd, was so intrigued with the basic 
concept that he pursued a thesis focused on the simple, 
reliable design of three satellite subsystems: structure, 

propulsion and attitude control. LCDR Boyd proposed a Shuttle 
GAS canister deployable satellite, like GLOMR or NUSAT. The 
GAS program is a Shuttle option that allows inexpensive, 
autonomous experiments to fly space available in the cargo 
bay [1], Due to the unceasing efforts of LCDR Boyd and Dr. 
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Fuhs, agencies such as DARPA, SDIO (Strategic Defense 
Initiative Organization) and Navy/Air Force STP (Space Test 
Program) Offices began showing a marked interest in the GAS 
satellite named ORION and in small spacecraft in general. 
This interest prompted the USAF STP office to pursue the 
development of an extended and enlarged GAS canister with an 
improved launch mechanism. The improved GAS can illustrated 
in Figures 2-5 and 2-6 can hold a satellite that is 35 inches 
long, 19 inches in diameter and can weigh up to 250 pounds 
[1]. In addition to the NPS satellite, three other efforts 
were launched in 1986 toward small, general purpose design. 
These included: "T-SAT" by INTRASPACE Company, an advanced 
NUSAT by GLOBESAT and American Rocket Company and an Orion 
type design by TRW [1]. This sudden flurry of activity in the 
small satellite business was capped by an encouraging DARPA 
program initiated in early 1988 called Lightsat. The Lightsat 
program promotes low-cost, lightweight satellites to be 
utilized by operational field commanders [1]. It appears that 
the idea of the small, general purpose satellite has finally 
come of age. 

C. ORION PRELIMINARY DESIGN 

The preliminary design of ORION adopted by LCDR Boyd is 
simple, reliable and relatively inexpensive with a price tag 
of about $1.5 million. ” A general purpose architecture was 
defined as the ability to successfully integrate various 
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payloads of the proper size while providing a propulsion, 
attitude control and standardized electrical, data and 
attitude control interfaces.” [1] 



GAS 

BEAM 




EXTENDED GAS 
CANISTER 



Figure 2-5 USAF/BALL Aerospace Extended GAS Canister [1] 
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General purpose naturally implies that trade-offs must be 
made in the design of the payload in order to reap the 
benefits of a low cost bus. Trade-offs are weight and size 
driven. The payload of any generic bus is constrained to 
weight and volume limits that do not affect most designs where 
the bus is built around the payload. A brief description of 
ORION subsystems will allow insight into the preliminary 
design . 

The ORION structural subsystem weighs approximately 40 lbm 
and consists of a .62" thick baseplate, three .75" honeycombed 
equipment decks, four extruded 1/16" thick longerons which 
house 78.5" long booms, a 1" thick honeycombed propellant tank 
strongback and four .05" thick semicircular skin panels. The 
baseplate, longerons, booms and skin panels are made of 7075- 
T6 aluminum while the honeycombed items use stainless steel 
facing material. The booms are designed to support 
magnetometers or 2 lbm tip weights to provide stability about 
the longitudinal spin axis. During Shuttle launch, this 
cylindrical structure will experience less than .035" of 
deflection. In compliance with NHB 1700. 7A safety 
specifications, the ORION model resonates above 35 Hz at modes 
of 160.8, 178.3 and 244.6 Hz [l]. Figure 2-7 demonstrates the 
basic structural subsystem. 

The propulsion subsystem uses hydrazine monopropellant and 
Shell 405 catalyst. The hydrazine tank is a spherical 
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Figure 2-7 ORION Structural Mock-up [1] 
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positive expulsion type. It is nitrogen pressurized and can 
hold up to 71.5 lbm of monopropellant . One 40 lbf thruster 
will be used for orbital maneuvers while six 0.1 lbf thrusters 
will be used for attitude control [1]. Design requirements 
included the ability to transfer from a Shuttle orbit of 135 
nm to 835 nm, the lower limit of the Van Allen radiation 
belts, in addition to a lifetime expectancy of 90 days to 3 
years [ 1 ] . 

The ORION power subsystem consists of 1765 in 2 of the 
outer skin covered by solar cells and a back-up supply of 24 
NiCad cells rated at 1.25 volts each [1]. An end-of-life DC 
requirement of 60 watts is based on the following assumptions: 
(1) one third of the skin surface area will be exposed to the 
sun, (2) solar cell efficiency of 14%, (3) solar constant of 
125.6 watts/ft 2 and, (4) 10% degradation per year. Even with 
these liberal assumptions, end-of-life power occurs after only 
18 months. 

The ORION attitude control subsystem consists of six 0.1 
lbf thrusters, one sun sensor, one horizon sensor and four 
magnetometers mounted on 78.5" booms. Spin stabilization 
about the longitudinal axis with ±1 degree of pointing 
accuracy was chosen as the simplest, most cost effective 
attitude control design option for ORION [1]. The fuel budget 
of the attitude control subsystem is limited by transfer 
altitude and effective deployment of the stabilizing booms in 
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this configuration. Figure 2-8 describes the attitude control 
and propulsion subsystems. 




Control Subsystems [1] 
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The remaining ORION subsystems were briefly discussed by 
LCDR Boyd. The thermal control subsystem consists of strip 
heaters and insulating blankets. The health and welfare of 
Orion will be monitored by a computer subsystem and telemetry 
link [1]. Note that the preliminary design assumes Space 
Shuttle deployment which is questionable due to the as yet 
approved use of hydrazine in a GAS can. After all the 
subsystem masses are tallied, a total of about 32 lbm remains 
for payload usage. Figure 2-9 shows the component placement 
for the ORION preliminary design. 
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Table 2-1 ORION COMPONENT MASSES [1] 



Component Mass(lbm) 



Top plate 


1.65 


Longerons ( 4 ) 


5.80 


Payload Midplate 


1.65 


Propellant Tank Strongback 


4.20 


Baseplate 


9.75 


Launch Restraint Pins (8) 


. 80 


Pressurant Mid-deck 


1.05 


Batteries/Boxes ( 4 ) 


28.0 


Earth Sensor 


.30 


Propellant Tank 


13.0 


Hydrazine Control Thrusters (6) 


9.18 


Orbital Transfer Thruster 


1.93 


Attitude Control and Payload Computer 


16.0 


Telemetry Transmitter/Receiver ( 2 ) 


10.0 


Pressurant Gas and Bottle (2) 


5.00 


Fill and Drain Valve (3) 


2 . 10 


Assorted Propulsion Tubing 


2 . 00 


Pyrotechnic Valve (2) 


1.00 


Power Conditioning Circuits (2) 


4.00 


Payload 


32.0 


Skin and Fasteners 


12.0 


Solar Cells 


.70 


Booms (4) 


4 .20 


Hydrazine Fuel 


71.5 


Magnetometers ( 4 ) 


2 . 00 


Total Mass 


245.81 
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III. SPECIAL THREE-AXIS CONSIDERATIONS 



A three-axis stabilized ORION will impose additional 
constraints on the structural, propulsion and power subsystem 
design compared to a spin stabilized version. A major 
consideration is the increased fuel usage for a mass expulsion 
attitude control system (ACS) . There are several options 
available to extend the ORION lifetime. Although not 
discussed in this thesis, but promoted as avenues of further 
research, there exist other energy efficient three-axis 
attitude control schemes. [2] 

A. ALTERNATE THREE-AXIS SCHEMES 

Momentum bias systems consist of a momentum wheel and 
either mass expulsion thrusters or magnetic torquers to dump 
excess momentum and control motion about the other two axes. 
The momentum wheel supplies gyroscopic stiffness which, in 
essence, cancels unwanted moments about one axis. Significant 
fuel savings result in the thruster moment bias design in 
comparison to the all-thruster system, but the added mass of 
the wheel must be compensated to maintain the same payload 
mass. Magnetic torquers use the torque produced when a 
magnetic dipole cuts through the Earth's magnetic field, but 
they are limited to low Earth orbits. They provide an endless 
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reservoir of control energy as long as current flows in the 
satellite magnets, but they are not preferred if fast reaction 
times are required. The necessity for propellent is 
diminished unless orbital transfer is desired. [2] 

Gravity gradient boom and thrusters or magnetic torquers 
afford a less expensive form of control since the boom is 
passive. Accuracy is reduced compared to the active momentum 
bias system, but fuel is still conserved. 

Reaction wheel systems operate independently about each 
axis just as a three-axis mass expulsion system does. [2] 
Thrusters are employed to dump excess wheel momentum so a fuel 
savings is realized when equated to the all-thruster system. 
The disadvantages of reaction wheel systems include the 
increased weight, complexity and cost of three wheels which 
could disqualify this system for a small satellite of ORION 
caliber. 

B. FUEL SAVING METHODS 

In addition to the minimum fuel-time control laws 
discussed and simulated in Chapters V and VI, fuel 
conservation can be practiced by restricting the employment 
of the 40 lbf orbital insertion thruster. Confining usage of 
the 40 lbf thruster to negation of the effects of aerodynamic 
drag could liberate most of the hydrazine for on-station 
attitude control. 
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The ability to transfer between circular orbits is a 
luxury for a small satellite, but AV required is lessened when 
the booster that carries ORION can insert close to the 
requested altitude. The original ORION designers assumed that 
the Space Shuttle would be the only practical launch option. 
In light of the stringent safety requirements for a man-rated 
payload, the shuttle seems to be the least feasible launch 
alternative when dealing with hydrazine monopropellant. Table 
3-1 lists current launch vehicles, payload weights, altitude 

and launch costs [6]. 

Table 3-1 SMALL SATELLITE BOOSTER OPTIONS 



Booster Name Payload (kg) Altitude (km) Cost(10 6 ) 



Titan 


14 , 030 


LEO 


$90+ 




4525 


GEO 




Atlas/Centaur 


5882-7625 


LEO ( i=28 . 5 ° 


$40 




2330-3145 


GTO ( i=28 .5° 




Delta II 


4460 


LEO 


$35-50 




1891 


GTO 




Scout G-l 


219 


555 


$11-12 


*Minuteman III 


680 


LEO (polar) 


$5 


*Conestoga II 


724 


LEO ( i=98 ° ) 


$7.03 




1040 


LEO ( i=38 ° ) 




*Super Starbird 


2080 


LEO ( i=0 0 ) 


$6+ 




724 


LEO (polar) 




Liberty 1A 


7260 


LEO 


$2-4 


Pegasus 


407 


LEO ( i=0 0 ) 


$6 




271 


LEO (polar) 




Foreiqn Systems 








Ariane IV(France 


4200 


GTO 


$80+ 


Long March (China) 


1267-2489 


LEO 


$35+ 




271-1358 


GTO 




*H-2 (Japan) 


8978 


LEO 


•? 




1990 


GEO 




Proton (USSR) 


20,000 


LEO ( i=52 0 ) 


$25-40 




2000 


GEO 





*not commercially available yet 
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If ORION must transfer to a higher orbit, the fuel budget 
will be driven by the AV required. The equation for the 
velocity of a circular orbit is 



V c = (^/R) % (3.1) 

where ji=GM e =3 . 986x10 s km 3 /s 2 , G is the Universal gravitational 
constant and M e is Earth's mass. The radius R=R e +h, R e =6378 km 
is Earth's radius and h is the satellite altitude. A Hohman 
transfer is the most efficient means of transitioning between 
coplanar orbits and is depicted in Figure 3-1. 




The total AV required for a Hohman transfer is 

A Vl = Vj - V el 
av 2 = V e2 - V 2 

AV t = AVi + AV 2 (3.2) 
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where AV X is the difference between the elliptical transfer 
orbit velocity at perigee and the inner circular orbit 
velocity and aV 2 is the difference between the transfer orbit 
velocity at apogee and the desired circular orbit velocity. 
The elliptical orbit transfer velocity is 

V = ( ( 2/i/R) + E) % (3.3) 

The specific mechanical energy of the transfer orbit is 

E = - p/ 2a (3.4) 

where a is the length of the semi-major axis of the transfer 
ellipse. The thrust equation for a rocket engine is 

F = m e V e + A e (p e - p 0 ) (3.5) 

m e = mass flow rate 
V e = exhaust velocity 
A e = exhaust cross-section 
p e = exhaust pressure 
p 0 = ambient pressure 

Specific impulse is a performance measure of the thrust to 
mass flow rate ratio for a rocket 

ISP = F/m (3.6) 
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ISP is a function of propellant type. The specific impulse 
of hydrazine is ISP=235 sec. The total impulse is the 
aggregate thrust obtained over time 



I T = J t Fdt = ISP Trip (3.7) 



where m p is propellant mass. The total propellant mass 
required to boost a satellite to a higher orbit using a Hohman 
transfer is 

m p = m 0 (1- e ( ’ AV/ISP8) ) (3.8) 



where m Q is the original satellite mass. The curves in Figure 
3-2 demonstrate the propellant mass required to boost ORION. 




Figure 3-2 Hydrazine Required by ORION for 
Transfer Orbits 
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ORION will only carry 71.4 lbm of hydrazine so the fuel 
remaining after transfer is alloted for attitude control. [1] 

C. SOLAR PANELS 

Three-axis stabilized satellite power subsystems generally 
include arrays of solar panels that extend after the transfer 
orbit. Panels allow more direct utilization of the sun's 
energy by maximizing the perpendicular rays of received 
sunlight through proper panel rotation. Solar cell output 
voltage is optimized in this configuration. This is the topic 
of another thesis design study and will not be examined in 
detail in this or later chapters. Solar panels will add more 
mass to the ORION structure and thereby add further size 
constraints to other subsystems. Suggested solar array 
designs for a three-axis ORION consist of cylindrical body 
mounted arrays that deploy or accordion style arrays that 
unfold and extend as shown in Figures 3-3 and 3-4. 




Mounted Solar Array [5] 
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Solar panels alter the moments of inertia of a satellite 
when deployed. ORION attitude control design assumes that 
I X =I Z which greatly simplifies the control problem by causing 
one axis to act independently of the other two in all regimes 
of motion. For I x *I y *I z , coupling occurs between all axes when 
the angular rate errors are large. An axially symmetric array 
would alleviate the coupling in one axis. 

The additional constraints on fuel and power budgets 
levied by a three-axis satellite are not insurmountable, but 
they must be investigated in more detail in future ORION 
design analyses. 
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IV. THREE-AXIS ATTITUDE DETERMINATION 



Satellite attitude control systems (ACS) consist of 
angular position and rate sensors or estimators and control 
actuators as shown in Figure 4-1. The sensors and estimators 
feed back angle and angle rate information in the form of 
error signals to be processed as required through a programmed 
control law. The actuators then react to commands from the 
ground or on-board microprocessors to implement corrections 
in attitude. The satellite is affected by unwanted internal 
and external torques such as gravity gradient, aerodynamic 
drag and thruster misalignment. An ACS conforms to the 
equations of motion and reacts to known and unknown 
disturbances. [2] 




Figure 4-1 Attitude Control Block Diagram [2] 
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A. SATELLITE DYNAMICS 

In modeling the ORION ACS, rigid-body dynamics will be 
assumed. Rigid body dynamics are used to approximate systems 
that have limited flexibility. The equations of motion for 
a satellite can be described by the movements of a point mass 
in space as shown in Figure 4-2. 




The mass P can experience rotation or translation with respect 
to the inertial reference frame XYZ. The position vector of 
P is 



R = Rc + r (4.1) 

where R c is a vector from the XYZ inertial frame to the xyz 
body frame and r is the position vector of P in the xyz frame. 

r= xi + yj + zk (4.2) 
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The absolute velocity of P is 



V = R = Rj. + r (4.3) 

where r includes linear and angular velocity in the body 
frame . 

r = r rel + w x r (4.4) 

Substituting (4.4) into (4.3), the absolute velocity becomes 

V = V c + V rel + u x r (4.5) 

The absolute acceleration of P is the derivative of absolute 
velocity. 



a= a c + a rel + 2« x V rel + w x (w x r) (4.6) 

In order to simplify future equations, two assumptions 
will be made. The xyz frame is internal to the spacecraft 
body so that v rel =o and the spacecraft center of mass (CM) 
coincides with the center of the body frame so that J m rdm=0. 
The linear momentum of a rigid-body is defined as the integral 
of the absolute velocity times the differential mass. 
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P = J~m(V c + w x r)dm = mV, 



(4.7) 



Angular momentum of a rigid-body is defined as the moment of 
the linear momentum or the cross product of r and P. 

H c = J m (r x V)dm = /„ r x (a, x r)dm (4.8) 

The angular velocity vector is equal to 

o> = w x i + w y j + w z k (4.9) 

Moments and products of inertia are analogous to the mass 
of a linear system but relate to mass distribution for a 
rotating rigid body. Inertia is the sum of the products of 
the component masses and their distance form the chosen axis 
of rotation. 

I = s mi r L z (4.10) 

Moments of inertia become 

I x = s iMy 2 + z 2 ) (4 . ll) 

I y = S m t (x 2 + z 2 ) 

i, = s %(* 2 + y 2 ) 
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Figure 4-3 shows ORION fuel and dry mass densities along the 
X and Y axes. 
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Table 4-1 lists preliminary subsystem masses and principal 
moments of inertia for ORION. 



Table 4-1 SUBSYSTEM MOMENTS OF INERTIA [12] 



Subsystem Mass (kg) lx Iy Iz 



Structure 


28.99 


2.4652 


2.4652 


1.4148 


Propulsion 


44.57 


1.6218 


1.6647 


0.8844 


Telemetry 


4 . 61 


0.6089 


0.5500 


0.0652 


Attitude Control 


4.41 


9.3621 


9.3827 


18.6291 


Computer 


3.45 


0.3482 


0.3482 


0.0214 


Electrical Power 


13 . 84 


0.7379 


0.5016 


0.3011 


Thermal 


0.0 


0.0000 


0.0000 


0 . 0000 


Payload 


11.22 


1.0812 


1.1117 


0.0916 


Total (kg-m 2 ) 


111.09 


16.23 


16.02 


21.04 


Total (Slug-ft 2 ) 


7.61 


11.9 


11.82 


15.97 


Without booms 










Total (kg-m 2 ) 


106.68 


6.8633 


6.6415 


2.7758 


Total (Slug-ft 2 ) 


7 .31 


5.06 


4.90 


2.05 



The products of inertia below are zero when the principal 
moments of inertia are chosen as rotation axes. 



Ixy = 2 m, xy 
I xz = E m A xz 
I yz = S m A yz 

After solving (4.8) and recognizing (4.9) and (4.11), the 
angular momentum becomes 





I* "~Ixy "”^xz 




’«x’ 


H c = 


i 

H 

a 

H 

1 

H 

N 




u y 




_ Ixz ""Iyz Iz _ 




_ W Z. 
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The inertia matrix is called the inertia tensor. To simplify 
the design of ORION, the inertia tensor will be assumed 
diagonal with zero cross terms. 

The kinetic energy of a rigid body is defined as the 
integral of the dot product of the absolute velocity. 

T = VV dm (4.13) 

This thesis will focus on rotation only. Substituting (4.5) 
into (4.13) the kinetic energy reduces to 



T= Z [ Wy W 2 ] 


[1,0 0 1 








O 

H 

O 




“y 




.0 0 1 z. 







Newton's Second Law can be applied to a rigid body as 

F = / m adm = ma c (4.15) 

A moment about the satellite CM is the cross product of 
distance and force. 



M c = H c + J m R c x rdm = H c 
= H crel + o> x H c 



(4.16) 



Angular momentum in the body frame can be expressed as 
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H c = H x i + Hyj + H,k 



(4.17) 



where xyz coincides with the principal moments of inertia. 
From (4.12), the elements of the angular momentum vector are 

H x = I x w x (4.18) 

Hy = Iy“>y 

Hz ^z^z 



After some manipulation, (4.16) is transformed into the Euler 
Moment equations. 



M x = 




+ 


w y w z ( iz 


- Iy) 


(4.19) 


M, = 


• 

J y w y 


+ 




" I z ) 




M z = 


1 z u z 


+ 


W X Wy(! y 


" I x ) 





Equation (4.19) will be used extensively through out this 
thesis to simulate ORION attitude control. 

An inertial coordinate system with X 0 Y 0 Z 0 fixed at the 
center of the Earth as shown in Figure 4-4 is used as a 
reference frame when describing satellite motion. The roll 
axis, X, describes the orbital velocity vector with magnitude, 
u> or which is assumed constant for a nearly circular orbit. 
The yaw axis, Z, extends from the satellite CM toward the 
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Earth CM. The pitch axis, Y, is perpendicular to the orbital 
or XY plane. Errors caused by disturbance torques are 
described by the angles <f > , V and 8 about the roll, yaw and 
pitch axes, respectively. In unperturbed form, the Euler 
moment equations become 

= (ix ~ (4.20) 

I y w y = (I 2 - I x )w x w 2 

I z“z = (lx - iy) u x w y 
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The right hand side of (4.20) exhibits the coupling that 
occurs when angular error rates are large. For on orbit 
control of a three-axis stabilized system the angular 
velocities are assumed small which reduces (4.20) to 



» 

3> x = 0 (4.21) 

IyWy = 0 

I 2 w 2 = 0 

Satellite acquisition is a stage in which initial large 
attitude adjustments are made to correct for booster 
separation. In this phase the satellite frequently 

experiences a tumbling motion which translates to large errors 
in angular position and rate. Excess rates must be negated 
before the on orbit control law can be implemented. Angular 
velocity errors as viewed from the inertial frame can be 
expressed as 



w x = <f> - w 0 V> (4.22) 

Wy = 6 ~ 

w 2 = V + w 0 <£ 

Replacing rate errors in the Euler moment equations yields 
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M x = I X U - u 0 i>) + (I, - I y ) w«>(tf + u 0 4>) (4.23) 

My = ly * 

M z = I,(tf + » 0 4 >) + (I y - I x )« 0 (* - co 0 ^) 

B. EXTERNAL DISTURBANCES 

The major external disturbances for low altitude orbits 
are caused by the Earth's gravitational field and aerodynamic 
drag. Solar pressure contributions are considerd negligible 
at low altitude but a dominant factor beyond geosynchronous 
altitudes . 

1. Gravity Gradient Torque 

Gravity gradient torques can be used as a passive form 
of three-axis stabilization for Earth pointing payloads if a 
boom is attached to the satellite. This torque results from 
the interaction between the earth-satellite masses. [2] The 
gravitational moment can be expressed as 

Mq = J ra r x F dm (4.24) 



where the gravitational force is 



F 



M (R e ~ r) dm 

R. 3 



1 + 



3r • Rg 
R } 



(4.25) 



Substituting (4.24) into (4.23) results in 
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Mg = 



(4.26) 



— f(r X Re) (r • RJ 
Re' 



where R e , the vector between Earth-satellite centers of mass 
as shown in Figure 4-5 is expressed in the satellite body 
frame as 



R e = Re(~ s i n ^i + sin^costfj + cos<£cos0k) (4.27) 
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For small angle corrections equation (4.26) reduces to 



2 . Aerodynamic Drag Torque 

Aerodynamic drag is caused by the friction between the 
satellite as it orbits at an angular rate of w 0 and the 
Earth’s atmosphere. It is felt most below 500 km. The 
effects of drag cause a general decrease in angular velocity 
thereby resulting in decreasing altitude or orbital decay. 
The aerodynamic drag force is a function of air density, p ; 
velocity, V; surface area perpendicular to the flight path, 
A; and a dimensionless number called the coefficient of drag, 
C D . [8] 



The orbital lifetime of a satellite is a function of altitude, 
h, and eccentricity, e. [7] 




(4.28) 



i pV 2 C d A 



(4.29) 



6378 + hi 3/2 



(4.30) 



L = NP = 1.15185X10' 7 N 



1-e 
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where N is the estimated number of revolutions and P is the 
period of the orbit. The values of N can be obtained directly 
from Figure 4-6 where the ballistic coefficient for ORION is 
m/ (C D A) =132 . 87 kg/m 2 . Note that the eccentricity for a 
circular orbit is zero, e=0. For ORION, orbital lifetime is 
shown in Figure 4-7 as a function of altitude. The linear 
velocity of a satellite equals 

V = w 0 R = (fi/R) k (4.31) 
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ORION LIFETIME 




Figure 4-7 ORION Lifetime due to 

Aerodynamic Drag 



where R = R e + h , the altitude above Earth CM. The drag 
coefficient is between 1 and 2 in SI units. The value used 
for ORION will be C D =2 . The surface area of ORION 
perpendicular to the flight path is A=.429m 2 . This reduces 
(4.29) to 



F d = 1.71 x 10 n p/R (4.32) 

ORION drag torque will can be summed up as 

Tdrag = F D L cosV- (4.33) 
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where L=.076m is the ORION center of pressure location along 
the y axis. [8] The effects of aerodynamic drag and 
gravitational gradient are simulated for a three-axis ORION 
in Chapter V. 

C. SENSORS AND ACTUATORS 

The ACS sensor suite is dictated by the type of 
stabilization, pointing accuracy, altitude, orbital parameters 
and cost constraints of the satellite budget. There are 
potentially three distinct stages of attitude control that 
ORION must handle. These include transfer orbit, acquisition 
and small error correction which includes altitude 
maintenance. If a transfer orbit is required, the satellite 
will most likely be spinning during this phase so different 
types of attitude determination sensors will be needed to 
accommodate spinning and non-spinning modes. In the 
acquisition stage, the satellite experiences large error rates 
which must be negated before accurate pointing can be 
sustained. The small error stage is the steady state 
operation that the satellite will undergo for the majority of 
its active life. 

The attitude of a spacecraft in inertial space can be 
completely realized by two vectors as long as they are not 
colinear. Point sources such as the sun and stars afford 
excellent reference vectors while the Earth and planets 
present extended sources. Pitch, roll and yaw attitude errors 
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relative to the desired satellite orientation must be 
determined in addition to error rates about these three axes. 
Combinations of horizon and sun sensors can furnish position 
error while rate sensing gyros or error estimation routines 
can be used to resolve rate errors. 

The current selection of attitude sensors for a spinning 
ORION consists of sun sensors, horizon sensors and 
magnetometers. The sun and horizon sensors provide the 
primary angular error output while the magnetometers act as 
back-up sensors. The spinning motion of the original 
satellite simplifies the number and cost of primary sensors 
by producing a periodic search motion. In a three-axis 
stabilized system this scanning motion must be integrated into 
the instrument in the case of the horizon sensor while 
multiple sun sensors with wider fields of view must be 
employed. 

1. Horizon Sensors 

Earth or horizon sensors are infrared devices that 
detect the 14-16 |im radiance profile of the earth as the space 
to earth and earth to space transitions occur. Horizon 
sensors can be used to ascertain pitch and roll errors. [2] 
Scanning or radiation balance detectors are employed in three- 
axis designs to provide a sweeping action across the disk of 
the Earth. Scanning sensors use two beams centered on a 
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reference line to scan the Northern and Southern hemispheres 
in an east-west direction. The scan pulse width is 

pw= d/(2hw h ) (4.34) 

where d is the distance between horizon crossings, h is the 
satellite altitude and w h is the scanning rate of the sensor. 
The difference between in-crossing to reference and reference 
to out-crossing pulse widths for each beam can be used to 
determine pitch error while the difference between Northern 
and Southern scans is proportional to roll error. An 
incremental pulse generator outputs a pulse at fractions of 
a degree for the known scan rate. Pulse width is computed by 
adding pulses from acquisition of signal (AOS) to reference 
then from reference to loss of signal (LOS). [7] Typical 
output from a scanning horizon sensor is depicted in Figure 
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2 . 



Sun Sensors 



Since two non-coincidental directions are needed to 
measure alignment uniquely, it is necessary to find another 
reference vector besides that from the Earth. The sun or 
another star are logical alternatives since they can be 
approximated as point sources. Stars are the best inertial 
reference points and are not coplanar unless the inclination 
is zero and the orbital path is in the equatorial plane. Star 
sensors are orders of magnitude more complex, expensive and 
bulky than sun sensors which owe their simplicity to the close 
proximity of the most radiant body in the solar system. 

Sun sensors are the most commonly used attitude sensor 
due to their relatively low cost, lightweight and broad 
application. [2] They can render accurate position error 
updates in harmony with horizon sensors when the satellite to 
Earth and satellite to sun vectors are non-parallel. Vector 
alignment can occur up to twice an orbit at local noon and 
midnight. Sun sensors are basically silicon solar cells with 
specialized masks for varying functions. Analog and digital 
sensors are available to measure the angle between the sun 
line and the pitch axis or sun presence in a particular field 
of view for instrument protection. The sun line vector is 
defined as 



S = sinacos5I +sin6J cosqcosSK (4.35) 
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where 8 is the sun declination and a is the orbital angle of 
the satellite measured form local noon. Yaw errors about the 
pitch axis in the inertial frame can be expressed using 

j = — V>I + J + <f>K (4.36) 

for small angles. The sun angle, 6 S , is obtained by taking 
the dot product of the sun vector and j 

cos0 s = S • j = -V»sinacos5 + sins + ^cosqcosS (4.37) 

As the yaw error changes the sun angle varies 

dCOS0 s 

= -sinacosS 

av> (4 ‘ 38) 

thus yaw error measurement is obtained. [2] 

3 . Magnetometers 

Magnetometers are not a required sensor for current 
or future ORION designs. The original design utilizes booms 
with 2 lbm weights at the end to enhance the stability of the 
spinning satellite by altering the principal moments of 
inertia. Three of the tip weights are magnetometers which 
will provide a back-up attitude determination capability that 
dead weights cannot. Magnetometers are solenoid coils that 
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obey Faraday's Laws of Magnetic Induction. When placed in a 
time-varying magnetic field, solenoid coils will produce a 
voltage 



V = f E ■ dl = -AN/x (dB A /dt) (4.39) 

where E is the electromotive force, A is the cross-sectional 
area of the solenoid core, N is the magnetic permeability of 
the core and Bj^ is the magnetic field component along the 
solenoid axis. The voltage is proportional to the magnetic 
field magnitude and direction. On-board control using 
magnetometers requires knowledge of the Earth's magnetic field 
at the orbit altitude. Knowledge of the magnetic field is not 
complete at low altitudes and too weak above 1000 km since the 
field strength decreases as the cube of the distance from the 
Earth's CM. [7] 

4 . Rate Gyroscopes 

In the acquisition phase, some means of determining 
error rates is required before they can be cancelled. Rate 
sensing gyroscopes, Luenberger observers or Kalman filters for 
the sun and horizon sensor output are viable solutions. Rate 
integrating gyros (RIG) in triplets are the most accurate 
instruments used to detect error rates and positions about 
three axes. The output of a RIG is the current required to 
maintain the gyro gimbals at a null or zero position. The 
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current is proportional to the torque needed to zero the 
rates. The total angular momentum of a gyro is 



H = LS + I o 0O (4.40) 

where L is the magnitude of the angular momentum of the rotor 
about the spin axis and I Q 6 is the magnitude of the angular 
momentum of the gimbal system about the output axis. The 
nulling torque is 



T g = (dH/dt) gyro + « x H (4.41) 

Angular velocity in the gyro inertial frame is 

u) = (jjl WqO w s S ( 4 . 42 ) 

where I is the gyro input axis, S is the rotor spin axis and 
O is the output axis. Figure 4-9 depicts gyroscope operation 
and the physical location of the axes described by equations 
(4.40) and (4.42). The differential equation of motion for 
a gyroscope is 

• • • 

I o g + D9 + KS - « X L = 0 (4.43) 

where D is a viscous damping term and K is a restoring term. 
[7] 
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5. Hydrazine Thrusters 

The functions of station keeping and attitude error 
correction are performed solely by control actuators called 
thrusters or reaction jets in this design. The thrusters on 
ORION use the monopropellant hydrazine with the Shell 405 
catalyst for these functions. The hydrazine storage tank is 
pressurized with nitrogen. A basic hydrazine reaction jet is 
pictured in Figure 4-10. 

The lifetime of a three-axis stabilized satellite is 
more acutely limited by the operational usefulness of the 
Shell 405 catalyst than in a spin stabilized vehicle where 
fuel usage can be more accurately predicted. The reaction 
jets on a three-axis stabilized satellite will be subjected 
to intermittant firing or low duty cycle firing. Thruster 
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Direction of thrust 




Nozzle 



Figure 4-10 Basic Thruster Design [2] 

performance is largely determined by the temperature of the 
catalyst bed. During steady state operation, intermittant 
firing of the jets will not allow the catalyst bed to maintain 
sufficient warmth to facilitate the most efficient thrust 
profile between firings. Typical values of specific impulse 
for hydrazine during transfer and steady state operation are 
235 and 135 seconds, respectively. The specific impulse is 
seriously degraded in the steady state regime. Ultimately, 
this leads to the conclusion that the catalyst bed in a low 
duty cycle environment must be heated which requires 
considerably more power than just that used for opening and 
closing valves and heating the hydrazine to prevent freezing. 
[ 2 ] 

In addition to the catalyst bed temperature, thruster 
performance is affected by impurities in the hydrazine that 
accumulate in the catayst bed after each firing and variation 
in the nitrogen pressure as the propellant is depleted. [2] 
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Figure 4-11 shows the linear relationship between thrust 
obtainable and pressure. 




A typical thruster profile is illustrated in figure 
4-12. This profile assumes that physical conditions are 
acceptable, but that there is a finite delay between the time 
that the thruster is commanded on and the time that it reaches 
maximum thrust and likewise for the command off sequence. The 
time from t x to t 2 and t< to t 5 is a function of the valve 
opening and closing action. In Chapter V, the thruster 
profile is modeled as a pulse function for simplicity. 

The current thruster selection for ORION includes a 
40 lbf jet for orbital insertion and six .1 lbf attitude 
control jets. Figure 4-13 lists design specifications for 
these thrusters. 
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Design Characteristics 



□ 


Propellant 


Hydrazine 


□ 


Catalyst 


Shell 405/lCH 202 


□ 


Thrust, Steady State (lbf) 


40—11 


□ 


Feed Pressure (psia) 


450— 100 


1 □ 


Chamber Pressure Ipsia) 


223 — 61 


□ 


Expansion Ratio 


20:1 


: □ 


Flow Rate (Ibm/sec) 


0 17—0.054 


□ 


Valve 


Wright Components Single Seat 


□ 


Valve Power 


50 Watts @ 28 vdc and 75 C F 


□ 


Weight (Ibm) 


1.95 




Engine 


1.35 




Valve 


0.60 


Demonstrated Performance 




□ 


Specific Impulse (Ibf-sec/lbm) 


235 — 203 


□ 


Total Impulse (Ibl-sec) 


4.350 


□ 


Total Pulses 


2,258 


□ 


Minimum Impulse Bit ( lbf- sec) 


0.71 @ 450 psia & 20 ms ON 


a 


Steady-State Firing (sec) 


y7 



Figure 4-13a Rocket Research 40 lbf Thruster - 



Model 107B [1] 
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Design Characteristics 

□ Propellant 

□ Catalyst . 

□ Thrust, Steady State (Ibl) 

□ Feed Pressure (psia) 

' □ Chamber Pressure (psia) 

□ Expansion Ratio 

□ Flow Rate (Ibm/sec) 

□ Valve 

□ Valve Power 

□ Weight (Ibm) 

Engine 

Valve 



3 SO 2 26 




Hydrazine 
Shell 405 
0.252—0.042 
.420—70 
370—60 
100:1 

0 . 001 — 0.0002 
Wright Components Dual Seat 
9 Watts Max. at 28 vdc and 45 ®F 
0.73 

0.28 

0.45 



Demonstrated Performance 

□ Specific Impulse (Ibt-sec/lbm) . . 227 — 206 

□ Total Impulse (Ibl-sec) 35.625 

□ Total Pulses 410,000 

□ Minimum Impulse Bit (Ibf sec) 0.001 © 150 psia & 8 ms ON 

□ Steady-State Firing (sec) 64,800 

Figure 4-13b Rocket Research .1 lbf Thruster 
Model 103C [1] 
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V. MINIMUM FUEL-TIME PROBLEM 



Three-axis stabilization is inherently more expensive than 
spin stabilization due to the increase in complexity of 
attitude sensors and the energy necessary to maintain the 
desired pointing accuracy. For Orion this translates to more 
fuel needed once on-station than the original spinning 
version. This limits the total AV available for orbit 
transfer and the active on-station control life of the 
satellite. Therefore, the effective operational lifetime of 
the payload could be dictated by fuel constraints before other 
considerations. In order to conserve fuel and provide maximum 
operational lifespan, a minimum control effort solution using 
Pontryagin's minimum principle is explored in this section. 

[9] 



A. ON-STATION CONTROL LAW 

For small angular rate and pointing errors, examination 
of one axis applies to all three axes. Thrusters are the sole 
control device for this design study and act independently in 
this case. The simplified differential equation for one axis 
is 

li = T c (5.1) 
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where T c is the thruster control torque and I is the inertia 
tensor with the products of inertia equal to zero. The 
normalized state equations for one axis are 



~Xj ( t ) 




'0 


1' 




’Xl(t)' 


1 

X • 

to 

✓ — s 

ft 

1 




0 


0 




X 2 (t)_ 



0] (5.2) 

u(t) 

1 



where x 1 is angular position, x 2 is angular velocity and u is 
the control torque provided by the thrusters. The minimum fuel 
performance index or cost function is 

J(u) = St (1 +A|u|)dt (5.3) 

where A is a weighting factor applied to the control. As A-+0, 
the performance index becomes the minimum time solution and 
as A-+<» , the performance index approaches the pure minimum fuel 
solution. Since any minimum effort problem is a trade off 
between elapsed time and consumed fuel both should be 
considered and are represented by (5.3). The performance 
measure must be minimized to obtain the optimal solution for 
given initial conditions and final states, x(t f )=0. Every 
value of the weighting factor provides a minimal trajectory. 
The selection of the weighting factor is determined by 
reasonable transient time and fuel expenditure. The 
normalized admissible controls for thruster operation are 
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|u(t) I < T c /I = N 



(5.4) 



where T c is the maximum torque available from the reaction 
jets. The cost function for thruster operation is called 
' Bang-Off-Bang ' since the satellite drifts before a final 
control torque is applied to bring it to the desired final 
attitude. The Hamiltonian for (5.3) is defined as 

H (x , u, p, t) = 1+ A | u | +p 1 x 2 +p 2 u (5.5) 

where p(t) are Lagrange multipliers. The costate equations are 



Pi(t) = - 8H/dx 1 = 0 (5.6) 

p 2 (t) = - 3H/3x 2 = -Pi(t) 

which have solutions 



Pi(t) = c x (5.7) 

P 2 (t) = -Cit + c 2 



The Hamiltonian must be minimized on an optimal trajectory, 
so the form of optimal control is 
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u(t) 



' N 
0 
-N 



p 2 < -A 
-A < p 2 < A 
p 2 > A 



(5.8) 



where switching occurs at | p 2 1 =A . The effect that any known 
or unknown torque has on the satellite can be best illustrated 
in the phase plane. This is a plot of angular position versus 
velocity as a function of applied torque. There are 
essentially two switching curves for the given minimum 
trajectory. These correspond to a curve that is dependent on 
the weighting factor, A, and the minimum time curve as shown 
in the phase plane representation in Figure 5-1. 



SWITCHING CURVES 




Angle 



Figure 5-1 Double Switching Curve Control 
Law 

If the initial conditions place the satellite in region 1, a 
negative torque is applied until t x when region 3 is entered 
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and control is turned off. Positive control is applied at t 2 , 
when the satelite reaches region 2, which is the switching 
curve for the minimum time solution. The reverse is true if 
the initial conditions place the satellite in region 2. At 
t 1 the control goes off, uftj^O, and (5.6) becomes 

P 2 ( t x ) = A = -C 1 t 1 + C 2 (5.9) 

The position continues to change as a function of velocity at 
turn-off. After t 2 , control is turned on, u(t 2 )=0, and the 
satellite ideally moves to the desired position. At t 2 (5.9) 
becomes 



P 2 (t 2 ) = "A = -Cita + c 2 (5.10) 

Subtracting (5.10) from (5.9) 



(t 2 - t x ) = 2A/c x 



(5.11) 



Assuming that velocity is constant at turn-off, x 2 =0 

x 2 (ti) = x 2 (t 2 ) = x 2 (t) 

Solving the Hamiltonian at the switching times and 
substituting (5.6) for p x yields 
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H(t x ) = H (t 2 ) = 0 
H(t x ) = 1 - C 1 X 2 (t 1 )=0 

c = L_ (5.12) 

X 2 (tj) 



Substituting (5.11) into (5.12) 



(t 2 - tj) = 2 Ax 2 ( t x ) 



(5.13) 



Solving the state equation, x x (t)=x 2 (t) 



Xi(t 2 ) - Xjfti) = x 2 (t x ) [t 2 -tj 



(5.14) 



From t 2 to t f the satellite follows a minimum time trajectory 
where 



Xi(t) 



i x/(t) 
|N| 



(5.15) 



Substituting 



(5.13) 

Xi(t) 



and (5.15) 

, ( 4 A + 1) 
2 I N | 



into (5.14) 



x 2 2 (t) = 0 



(5.16) 



which is the control law for this minimum fuel problem. 
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Equation (5.16) becomes usable in all quadrants when 
x 2 2 =x 2 |x 2 | which takes the sign of x 2 into account. Figure 5- 
2 demonstrates minimum trajectories for various values of A. 
The initial conditions x x (0) =1 , x 2 ( 0) =0 were used to test for 
minimization of the selected cost function. As shown in 
figure 5-3, the minimum value of the chosen cost function 
occurs at A=2.25 when this value is substituted into equation 
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Figure 5-2 Parabolic Switching Curves for 
A=0, 1, 10, 100 



Figure 5-4a is a time response of angular position and rate. 
Figure 5-4b is the time response of control torque. Notice 
that thruster chatter occurs at the end of the cycle, since 
fine control near the desired attitude is not implemented. 
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WEIGIITINC VARIABLE vs COST 




A=( 1 + 



Figure 5-3 Minimum Fuel-Time Cost Function 
J - /(I + 2 . 25 | u | ) dt 





Figure 5-4 Time Response of States and 
Control 
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Figure 5-4b shows the "bang-deadzone-bang" behavior of the 
control law. Figure 5-5 is a phase plane representation of 
the same cycle. 




Figure 5-5 Phase Plane for A=2.25 



B. LIMIT CYCLE CONTROL 

Ideally, the last thrust period will take the satellite 
to the origin, but this is a function of the integration step 
size. A similiar analogy exists for the thruster in the form 
of the minimum impulse bit available. Physically angular 
rates are never zeroed, but they can be maintained in a 
region. Satellite motion is then dictated by a finite limit 
cycle even in the absence of external disturbances. 

While in the limit cycle attitude errors are reduced below 
specified values. The fine tuning region is determined by 
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pointing accuracy requirements, but it is also a function of 
minimum jet impulse for angular velocity and the precision or 
noise contribution of the attitude sensor instrumentation. 
Minimum jet impulse is a measure of the smallest thrust unit 
available when the control valves are immediately opened and 
closed, i.e., the shortest single jet pulse. In order to 
increase the coasting period in the limit deadzone, the 
satellite rate must be reduced to the lowest value possible 
without causing thruster chatter. Chatter is a constant 
cycling of opposing thrusters when limits are exceeded after 
each firing. The limit cycle for the cost function simulated 
in this example is depicted in Figure 5-6. 
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A third switching curve was implemented in the boundary 
region only to minimize fuel expended while cycling. The 
third switching curve has a flatter slope, A=250, and thus a 
wider deadzone. The best deadzone obtainable is again limited 
by the minimum jet impulse available from the thrusters. 
Figure 5-7 utilizes the third switching curve and can be 
compared to Figure 5-6 where only two switching curves were 
realized to demonstrate the improved fuel consumption. 




Figure 5-8 shows the time responses of thruster operation 
during fine tuning in both cases. Note the single impulse 
firing that occurs at the edge of the limit boundaries when 
the third switching curve is implemented. The appendix 
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provides a description of the Tutsim program that simulates 
the independent control law about a single axis. 







44 4G 40 50 

Time (sec) 
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Figure 5-8 Fine Control 



C. ACQUISITION PHASE CONTROL 

The acquisition phase is the most critical aspect of 
satellite operation after a successful launch. During 
acquisition, the satellite may be tumbling. Tumbling 
translates to high error rates which will cause coupling about 
at least two axes. In this phase, a minimum time solution is 
desired to cancel the excess rates quickly. [6] The Euler 
moment equations of motion are 
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(5.17) 



I X W X = (Iy - I Z )U y U> z + T x 

!y w y * (I* “ Ix)Wx W 2 + T y 

I z^z = Ux " I y) w x w y + T z 

where T x , T y and T z are the control torques; the product of 
thrust and moment arm. If a single axis of symmetry is 
assumed, as is the case for a cylinder, two of the moments of 
inertia will be equal 

lx = Iz = I 



The moment equations become 



w x = Qw y w z + u x (5.18) 

U)y — U y 

w z = -a w y w x + U z 



where a = (I-I y )/I and u^Ti/Ii. 

Motion about the axis of symmetry is independent of the 
other two axes and is controlled at all times as discussed in 
the previous section. Assuming that w y is constant when Uy is 
zero during rate reduction, the state equations become 



' X1 (t)' 




'0 ft 




’ x l(t)' 


. X 2(t)_ 




-ft 0 




_ X 2(t). 



1 0] (5.19) 

u(t) 

0 1 
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where B=aw 



x 2 =w z , 



The minimum time 



y, Xj— W X , 

cost function is 



u 1 =u x and u 2 =u z . 



J = ft dt 



(5.20) 



The Hamiltonian for the coupled system using both controllers, 
u x and u 2 , is 



H(x,u,p,t) = 1 + p x x 2 - p 2 x x + PiUi + p 2 u 2 (5.21) 

for 3=1. The costate equations are 

Pi(t) = - SH/dXj = p 2 (t) (5.22) 

p 2 (t) = - aH/ax 2 = -Pi(t) 

which have solutions 

Pi(t) = c 3 cos (t+5 ) (5.23) 

p 2 (t) = c 4 sin (t+5 ) 

To minimize the Hamiltonian on an optimal trajectory, the form 
of optimal control is 

Ui(t) = -sign (p 3 ) (5.24) 

u 2 (t) = -sign (p 2 ) 
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The control torque of each set of thrusters will reverse every 
half cycle while the control mix will change every quarter 
cycle. The minimum time switching curves can be derived using 



dx x x 2 + u x 
dx 2 -x x + u 2 

/(-Xi+UaJdXj = f (x 2 +u 1 )dx 2 

2 _ 2 

- u 2 x x _ - ^ 2 _ - Ul x 2 

2 2 



(5.25) 



Adding u^+Uj 2 to both sides of (5.25) and multiplying by two, 
yields 



(x x - u 2 ) 2 + (x 2 + Ui) 2 = u x 2 + u 2 2 (5.26) 

For normalized control torque, u 1 =±l, u 2 =±l, (5.26) becomes 

(Xi T l) 2 + (x 2 ± l) 2 = 2 (5.27) 

Equation (5.26) describes four circles with their centers at 
(±1,±1) in each quadrant. The switching curves are defined 
as the quarter of each circle that intersects the abscissa and 
ordinate in each quadrant as shown in Figure 5-9. Beyond 
these scalloped sections the switching curves are mirror 
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images of the half sinusoidal sections along each axis, but 
minimal error is experienced if switching is implemented about 
the axes themselves. 




Figure 5-9 Double Controller 

Minimum Time Switching Curves 



A phase plane representation is used to illustrate rate 
reduction where w x is the abscissa and w z is the ordinate. 
Figure 5-10 depicts the implementation of the four quarter 
section switching curves in a circular region with a radius 
of two for initial rate conditions of (±1.5, ±1.5). Figure 5- 
11 portrays switching about the coordinate axes in the region 
beyond a radius two. In the outer region, u x =-sign (w x ) and 
u z =-sign(w z ) . Figure 5-12 shows the control torque time 
response for u x and u z . 
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DOUBLE CONTROLLER 




-2 



-3 1 1 ' > ' 

-3-2-10 l 2 3 

Figure 5-10 Inner Region Switching 




Figure 5-11 Outer Region Switching 

a ) ct) x =(x) z =6 b ) cj x =— 7 o) 2 =5 
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Figure 5-12 Outer Region Control 



Since the system is coupled, only one controller is needed 
to reduce the rates about both axes. If u 1 =0, equation (5.25) 
becomes 

2 2 

x i - U 2 X X = x 2 (5.28) 

2 2 



Adding u 2 2 to both sides and multiplying by two 

(x x - u 2 ) 2 + x 2 2 = u 2 2 (5.29) 

which becomes 



(x x I l) 2 + x 2 2 = 1 



(5.30) 
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for u 2 =±l. Conversely, for u 2 =0 equation (5.30) becomes 



Xl 2 + (x 2 2 ± 1) = 1 (5.31) 

for U!=±l. Equations (5.30) and (5.31) describe two symmetric 
circles with centers at Xl =±l and Xz =±l, respectively. The 
switching curves for the minimum time solution are the half 
circles shown in Figure 5-13. Figure 5-14 demonstrates rate 
reduction for a single controller system. For initial rates 
in the outer region, switching about the axes was implemented. 
Control torque values of u z =sign(w x ) for (5.30) and u x =sign(w 2 ) 
for (5.31) were applied. Figure 5.15 demonstrates large rate 
error reduction for both single controller cases. 




Curves 
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ffx 




Wx 

Figure 5-14 Single 




Figure 5-15 Outer 
Region w x =w z =6 
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Figure 5-16 illustrates the control torque time responses in 
the outer region. Comparison of Figures 5-12 with 5-16 shows 
that the double controller system takes half the time using 
the same amount of control effort as the single controller 
system. 




0123450789 10 



Time (sec) 




Time (sec) 

Figure 5-16 Outer Region Control Torque 

D. LIMIT CYCLE CONTROL FOR ORION 

ORION specifications can easily be implemented in the 
simulation described in previous sections. A three-axis 
stabilized ORION will be oriented such that the yaw axis is 
col inear with the longitudinal axis. The symmetry of the 
principal moments of inertia becomes I x =I y vice I X =I Z as in the 
spin stabilized configuration. ORION principal moments of 
inertia from Table 4-1 are 
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I x * I, « 6.753 kg-m 2 
I z = 2.776 kg -in 2 

Attitude control thruster data is obtained from Figure 4- 
14b. The thrust available from the attitude control reaction 
jets is between .24 and 1.068 Newtons. The minimum impulse 
bit is .0223 N-sec for .025 seconds. The total pulses 
available as a function of the catalyst effective lifetime are 
161,000. The mass flow rate of the thrusters is between .104 
and .453 grams/sec. 

Torque is the product of thrust and moment arm. The 
length and radius of the ORION cylinder are .889 and .2414 
meters, respectively. Thrusters operate in pairs to prevent 
translational movement so the average torque produced for 
ORION is 



T x = T y = 2 ( . 445) ( .889) = .7912 N-m 
T z = 2 ( .445) ( .2414) = .2148 N-m 

The control accelerations are 

u x = u y = .1172 s’ 2 
u z = .0774 s' 2 
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Pointing accuracies of ±1° and ±.05°/sec were used in the 
example illustrated to define the limit region. The minimum 
fuel-time switching curves were implemented with the values 
of A previously used, A=2.25 and 250 for the outer region 
deadzone and limit region deadzone switching curves, 
respectively. Figures 5—17 and 5—18 illustrate the minimum 
fuel-time control law for yaw and roll-pitch when ORION values 
for torque and moments of inertia are used. 




Figure 5-19 shows the yaw thruster operation for initial 
conditions of angle=.25 rad and zero rate. 
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Figure 5-18 Roll-Pitch Limit Cycle for 
Limits of ±1°, ±.05°/ sec 
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It is difficult to predict fuel usage up to the 
operational phase of flight for three-axis control. The 
approximate limit cycle lifetime with all thrusters operating 
is 



Total Fuel Available 

L = At (5.32) 

Fuel per Control Cycle 

where At = one limit cycle. [1] In the limit cycle, with no 
environmental disturbances added, individual pairs of 
thrusters pulse on four times per cycle. Using a mass flow 
rate of .104 grams/s for .025 sec, the minimum pulse time, 
gives a total fuel per control cycle of between .0624 and 
.0832 grams if all six thrusters are operating to maintain 
attitude in their respective limit regions. The yaw and roll- 
pitch control cycle duration for this example are 108 and 72 
seconds, respectively. Assuming that 2/3 of the original 
32.35 kg of hydrazine is available for limit cycle attitude 
control, total limit cycle lifetime would be about 320 days. 
Approximate limit cycle lifetimes for various combinations of 
weighting factors, rate limits and minimum impulse bits based 
on the hydrazine budget are listed in Table 5-1. Increasing 
the value of A or decreasing the rate error limit will 
decrease the rate error and thus lengthen the thruster off 
time, but the maximum effective value of A is limited by the 
the minimum impulse bit on-period. For the thruster used 
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Table 5-1 LIMIT CYCLE LIFETIMES 



u z At z (sec) rate(°/sec) A At cz (sec) L(days) 



.0774 


. 025 


. 075 


100 


48 


143 


.0774 


. 01 


.05 


250 


108 


323 


. 0774 


. 01 


.025 


1000 


125 


427 


.0774 


. 006 


.015 


1000 


185 


632 


. 0387 


. 025 


.075 


100 


86 


257 


. 0387 


. 025 


. 075 


1000 


94 


281 


. 0387 


. 01 


. 015 


1000 


218 


869 


Uxy At xy (sec) 


rate ( °/sec) 


A At c3cy (sec) 




. 1172 


. 025 


. 085 


100 


30 




. 1172 


. 025 


.05 


250 


72 




. 1172 


. 007 


.025 


1000 


112 




. 1172 


. 004 


.015 


1000 


180 




. 0586 


. 025 


.075 


100 


58 




.0586 


. 025 


.075 


1000 


64 




. 0586 


. 008 


.015 


1000 


195 





in this example, A>100 has negligible effect since the minimum 
AV imparted is 

AV Z = u z At = .0774 (.025) = .001935 rad/sec 
AV^ = u^y it = .1172 (.025) = .00293 rad/sec 

which will kick the rate error from the negative limit of 
-8.73x10’* rad/sec to a positive value of .0012 rad/sec. The 
example illustrated in the figures uses an integration step 
size of .01 seconds. In order to implement the fuel saving 
choices listed in table 5-1 attitude control thrusters with 
smaller minimum impulse bits and\or on-times must be utilized. 
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Limit cycle lifetime based on the number of available 
pulses ranges from 30 to 90 days. This is considerably less 
than limit cycle lifetimes based on fuel consumption. In 
addition to smaller impulse bits, the number of available 
pulses must increase to lengthen the active life of ORION. 
Another option would be to operate for specified periods 
during each orbit vice continuously. Rather than replacing 
the thrusters with more complex and expensive items, this 
could increase the active lifetime of the existing design to 
acceptable levels. 

E. ENVIRONMENTAL DISTURBANCE EFFECTS 

Since the majority of the satellite lifetime will be 
spent in the operational or steady state mode, this is the 
phase when the effects of the external disturbances of gravity 
gradient and aerodynamic drag will play the largest part. 
Substituting the small angle approximation for the angular 
velocity vector and adding a disturbance moment to the Euler 
Moment Equations yields 



4> ~ a u a Z 4> + U x + D x (5.33) 

e = u y + D y 
0 = U 2 + D z 
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where a = (I y -I z )/I x , D A = T^/Ii is the disturbance contribution 
and u A = Tci/Ii is the control offering. The small angle 
approximations for the gravity gradient moment equations are 



T g 



3w r 



f *<I,-Iy)l 

l o J 



(5.34) 



where the orbital rate w 0 2 = n e /R 3 . Figure 5-20 shows the 
orbital rate versus altitude for any satellite orbiting the 
Earth. 




Figure 5-20 Orbital Rate, w D , versus 
altitude 



Substituting (5.34) into (5.33) and eliminating the control 
term gives 
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4 > 



(5.35) 



<t> = "2u> 0 2 




e 





e 



which are equations of motion for simple harmonic oscillators. 
The moment of inertia ratios for ORION equal +.5929 and 
-.5929, respectively. Figures 5-21 and 5-22 illustrate the 
harmonic motion about the roll and pitch axes. Plots in 5- 
21a and 5-22a are phase plane representations while graphics 
in 5-21b and 5-22b show the time response for angle and rate. 
The motion is stable and periodic with a period that 
corresponds to the satellite orbit. The orbital period is 
defined as [4] 



Figure 5-23 furnishes the orbital period versus altitude for 
an Earth orbiting satellite. In Figures 5-21 and 5-22 the 
maximum rate error induced never exceeds the initial value, 
but the angle error will exceed the desired pointing accuracy 
over the course of an orbit if a rate error exists. Gravity 
gradient torque has a time constant that is equivalent to an 
orbital period. The time constant of this all-thruster 



P 




(5.36) 
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Figure 5-21a Gravitational Roll 
Effects Phase Plane 
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Figure 5-21b Gravitational Roll Effects 
Time Response 
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Figure 5-22a Gravitational Pitch 
Effects Phase Plane 





Figure 5-22b Gravitational Pitch Effects 
Time Response 
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EARTH SATELLITE ORBITAL PERIOD 




Figure 5-23 Orbital Period versus Altitude 



attitude control system is on the order of 100 seconds. If 
the control system is operating continuously to maintain the 
desired pointing accuracy, then the effects of gravity are 
insignificant compared to the normal cycling that occurs in 
the limit region. If there is zero rate error and the angle 
error is in the limit region, the gravity gradient torque will 
aid in maintaining the attitude in the limit region. 

The disturbance torques contributed by aerodynamic drag 

are 



T 



Drag 



-F d L sinflcostf’ 
0 

F d L COS8COSiJ> 



(5.37) 
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where F D is the drag force as a function of atmospheric 
density. Table 5-2 lists atmospheric densities at some 
altitudes. [8] 

Table 5-2 ATMOSPHERIC DENSITIES 

Altitude (km) Density (kg/m 3 ) 



100 


6xl0' 7 


200 


7xl0' 7 


400 


lxlO' 11 


1000 


lxlO’ 15 



At h=4 00 3cm the drag force is 2.52x10 A Newtons. For small 
angles (5.37) reduces to 



Dx — ~F d 


L 6 


= -1 . 92x10 5 


6 N-m 


Dz _ F d 


L 


= 1 . 92xl0’ 5 


N-m 



where L=.076m is the center of pressure offset from the center 
of mass. The effects of aerodynamic drag alone for an orbital 
period are illustrated in Figures 5-24 and 5-25 about the roll 
and yaw axes for initial conditions of angle=.01 rad and 
rate=.0001 rad/sec. The pitch angle error was assumed 
constant at .01 rad. The satellite response to drag is 
unstable, yet like gravity gradient, drag causes a slow 
response compared to the continuously operating all-thruster 
control system. Drag torques that would cause the attitude 
of the satellite to be unstable over the course of an orbital 
period do not affect motion during continuous control. 
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xIO-s* AERODYNAMIC DRAG RO!X EFFECTS 




Figure 5-24a Drag Roll Effects Phase Plane 
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Figure 5-24b Drag Roll Effects Time 
Response 
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AERODYNAMIC DRAG YA1T EFFECTS 




Figure 5-25a Drag Yaw Effects Phase Plane 
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Figure 5-25b 



Drag Yaw Effects Time Response 
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Finally the effects of gravity gradient and aerodynamic 
drag were combined. Substituting the disturbance torques into 
(5.33) yields 



4> = -a a> 0 2 <f> + U x - 2.84X10' 6 6 (5.38) 

which is the equation of a marginally stable oscillator. The 
effect that gravity and drag have on roll is illustrated in 
Figures 5-26 and 5-27 for initial conditions of angle=.01 rad, 
zero rate and rate=.0001 rad/sec, respectively. The 
oscillation does not increase without bound and actually never 
repeats the same cycle. The phase plane plots have a run time 
of 150 minutes while the time response plots have a run time 
of 300 minutes. Again, the time constant of the gravity-drag 
torque is orders of magnitude bigger than that of the all- 
thruster control system. If the control system were only 
operated at discrete intervals during each orbit the external 
disturbance torques could together run the rates up to 
unacceptable levels. 
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Zero Rate Error Phase Plane 





Figure 5-26b Gravity-Drag Roll Effects with 
Zero Rate Error Time Response 
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CRAVTTY-DEAC ROLL EFFECTS 




Figure 5-27a Gravity-Drag Roll Effects with 
Rate Error Phase Plane 




Rate Error Time Response 
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VI. SUMMARY AND CONCLUSIONS 



The purpose of this thesis was to design an all-thruster 
three-axis attitude control system for ORION, the Naval 
Postgraduate School satellite bus. Pontryagin's Minimization 
Principle was employed in the derivation of fuel-time 
minimization control laws for the acquisition and on-station 
phases of flight. The original design was modified slightly 
to accommodate a three-axis strategy. 

The original design history and specifications are 
discussed in Chapter II. The total mass of the ORION 
satellite is 250 lbm including 71.5 lbm of the monopropellant 
hydrazine for maneuvering. The attitude control system 
consists primarily of one 5 lbf thruster for orbital transfer 
and six .1 lbf thrusters for rotational control. The original 
mass placement and principal moments of inertia were left 
unchanged for simulation purposes. The four stabilizing booms 
are not utilized in this analysis. 

Additional concerns that affect a three-axis stabilized 
satellite compared to the original spin stabilized design are 
addressed in Chapter III. To meet power requirements, solar 
panels will be required vice just body-mounted solar cells. 
A three-axis ORION will be less fuel efficient than a spinning 
version. Various schemes using different combinations of 
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actuators are presented for future study. Alternate booster 
selections that could transfer ORION to higher orbits than the 
space shuttle and minimize the use of the 40 lbf thruster are 
suggested. 

The theory behind satellite motion, attitude 
determination and control hardware for a three-axis stabilized 
spacecraft, is developed in Chapter IV. Rigid-body spacecraft 
dynamics and the effects of the environmental disturbances of 
gravity gradient and aerodynamic drag are covered. Attitude 
control sensors and actuators are described with an emphasis 
on the higher degree of complexity required of these 
instruments vice those used on a spinning satellite. Sensors 
include horizon sensors, sun sensors, magnetometers, and rate 
gyros while the actuators are strictly thrusters. 

The derivation and simulation of the control laws that 
describe the acquisition and on-station phases of flight is 
found in Chapter V. It was determined that a dual-axis 
minimum time control law was optimal for reducing the rates 
as quickly as possible in the acquisition phase. Once the 
rates are reduced to acceptable levels, a fuel-time control 
law is implemented that utilizes a weighting factor to 
minimize fuel usage. 

Two regions of interest developed during the analysis of 
the fuel-time control law. The boundary region is a 
rectangular area centered about the desired position. The 
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deadzone is determined by the nominal pointing accuracy (±1° 
ORION design goal) in angle and the minimum thruster impulse 
bit for rate. The boundary region is distinguished by a four 
pulse limit cycle that repeats itself every 50-200 seconds for 
ORION. Assuming that 2/3 of the total hydrazine budget is 
available for limit cycle maintenance, lifetimes of 4-28 
months were obtained with all thrusters operating. The 
limiting factor for the selected thrusters was the number of 
pulses available before the Shell 405 catalyst lost 
effectiveness. The limit cycle lifetime decreased to 30-90 
days. The second region of interest was outside the boundary 
region, but close enough for the satellite to still experience 
independent or uncoupled motion about one or more axes. A 
double switching curve control law distinguished by "bang-off- 
bang" control was realized successfully. 

The environmental disturbances of gravity gradient and 
aerodynamic drag were also simulated in Chapter V. The 
gravity gradient disturbance imitates a simple harmonic 
oscillator with a period corresponding to an orbital period. 
A typical orbital period for a low earth orbit is on the order 
of 90 minutes. Aerodynamic drag applies an unstabe torque, 
but when coupled with gravity causes a random but marginally 
stable type of oscillation. Both disturbance torques are so 
small over the course of a limit control cycle, that the 
effects are negligible. 
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The three-axis stabilized design simulated in this thesis 
could be a viable attitude control option for the general 
purpose ORION bus. Upgraded thrusters must be selected that 
have smaller minimum impulse bit on-times and thrust values 
in addition to an increased catalyst lifetime if ORION is to 
operate in the constant update control mode for a feasible 
duration. Another possibility that is only mentioned in this 
work, concerns intermittent operation. ORION could function 
for specified periods every orbit when the desired target is 
in the field of view. This method would insure an extended 
lifetime, but the influence of external torques could be a 
limiting factor. 
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APPENDIX 



ATTITUDE CONTROL AND EXTERNAL DISTURBANCE PROGRAMS 



The following programs simulate a three-axis stabilized 
satellite in the acquisition and on-station phases. The 
effects of gravity gradient and aerodynamic drag are modeled 
in separate routines. 

The software package TUTSIM™ was used to execute the 
control laws derived in Chapter V [11]. TUTSIM™ models 
continuous time control systems in a signal flow format. The 
following commands are defined for the potential user: 



INT - integrator block 

CON - constant block 

GAI - gain block 

MUL - multiplier block 

DIV - divide block 

ABS - absolute value block 

AND - logical and block 

OR - logical or block 

IFE - If-then statement 

REL - relay block 



The numbers preceding commands are simply block identifiers. 
The numbers following commands are block inputs. The constant 
block does not have an input. The integrator and gain blocks 
require at least one input. The multiply, divide and logical 
operator blocks need at least two inputs. The if-then 
statement has three inputs. The first input is the argument. 
The second and third inputs are the commands that support true 
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or false arguments, respectively. The relay block has four 
inputs. The last input is the argument. The first three 
inputs are commands that support greater than, equal to or 
less than results, respectively. Further information on the 
capabilities of TUTSIM™ can be found in the operator manual. 

The first program simulates on-station control about one 
axis. Input parameters include initial conditions, principal 
moments of inertia, control torques, limit region boundary 
values and weighting factors for two switching curves. The 
command logic applies control according to the satellite 
attitude in the phase plane in relation to the switching 
curves and boundary limits. The control law for this program 
is derived in the first section of Chapter V. 

The next three programs simulate the acquisition phase 
for single and double axis control. Input parameters include 
initial conditions for rate errors, principal moments of 
inertia and torque values. These programs apply control as 
a function of the rates for coupled equations of motion. 
Switching occurs about the ordinate and abscissa in a defined 
outer region. Switching occurs on minimum time curves in the 
inner region. The control law for the acquisition routines 
can be found in the second section of Chapter V. 

The last two programs simulate the effects of gravity 
gradient and aerodynamic drag distubances, respectively. 
Inputs parmeters are provided for initial conditions, 
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principal moments of inertia and orbital altitude. The last 
section of Chapter V illustrates the effects of external 
disturbances . 
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